How the effective boson-boson interaction works in Bose-Fermi mixtures in periodic 

geometries 
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We study mixtures of spinless bosons and not spin-polarized fermions loaded in two dimensional 
optical lattices. We approach the problem of the ground state stability within the framework of 
the linear response theory; by the mean of an iterative procedure, we are able to obtain a relation 
for the dependence of boson-boson effective interaction on the absolute temperature of the sample. 
Proceeding from such a formula, we write down analyitical expressions for Supersolid (SS) and 
Phase Separation (PS) transition temperatures, and plot the phase diagrams. 



I. INTRODUCTION 

The possibility to achieve very low temperatures and 
the feasibility in the laboratory of implementing optical 
lattices represent the proper arena where testing the 
validity of certain condensed-matter theories and ob- 
serving the manifestations of effects related to quantum 
mechanical statistics. Within the interesting sinergy 
between quantum atom optics and condensed matter 
physics, recent trends are study of the superfluid to 
Mott-insulator phase transition in bosonic systems, (see 
EH and i), and the striving for the realization of a 
BCS-type condensate in a fermionic system (@ and 0]). 
When atomic systems made up of bosons and fermions 
are considered, a very rich scenario opens up (0], [H, 
*5|). In particular, Bose-Fermi mixtures sympathetically 
cooled into their quantum degenerate states (Q and 
M) exhibit a strong tendency towards Phase Separation 
( [1(1 and [ll|), demixingof the bosons and the fermions, 
and Supersolid ([III, E3, [HI)- Recently, these issues 
were addressed in |l2| . where spinless bosons and 
spin-polarized fermions confined in two dimensional 
(2D) periodic geometries were taken into account. As 
explained in [12j |. the effective interaction between the 
bosons of the mixture plays a crucial role to the end of 
determining the boundaries between the various phases 
of the system. 

The dynamics underlying Phase Separation is driven 
by a small change (Sns) in the bosonic density, which 
produces a modulation (Snp) of the fermionic density, re- 
lated to Sub by 5uf ~ — UBFN(Ep)5riB, with Ubf the 
on-site boson-fermion interaction amplitude, and N(Ep) 
the density of states at the Fermi energy Ep As 
consequence of the feedback of the fermionic distortion 
Snp, a shift of the bosonic energy — U B pN(Ep)(5nB) 2 /2 
occurs, thereby inducing an attraction between the 
bosons with strength U B pN(Ep). Phase Separation 
emerges when the induced attraction and the intrinsic 
repulsion between the bosons have the same order of 
magnitude [nl ]. 

One of the peculiarities of optical lattices is their 
periodic arrangement, that is the same of the crystal 
structure of a solid. The phenomelogical interpretation 



of particle localization at fixed sites at the base of Mott- 
insulator mechanism is related just to the crystalline 
structure . In Supersolid, such a structure is combined 
with the essence of a superfluid, i.e. stiffnes allowing for 
particles current to flow without dissipation. 
In general to supersolids are associated two kinds of 
order, which usually appear in mutual competition. 
These are the diagonal long-range order (DLRO) 
associated with the periodic density modulation in a 
crystal, and the off-diagonal long-range order (ODLRO) 
associated with the phase order in the condensate [20| . 
In the strongly interacting case, as in 4 He systems, 
(dH and [Hj]), supersolids have been proposed to exist 
and analyzed numerically in various model systems 
describing interacting bosons on a lattice [23[ . 

In this paper, we consider mixtures of spinless 
bosons and not spin-polarized fermions, so that a 
non zero s-wave interaction between fermions on the 
same site emerges. The goal of the present work is 
analyzing the effects of such an interaction on phase 
diagrams of the system. Approaching the pro blem 
within the linear response theory framework fl6l.ll8l [l9| , 
we show that it is possible to calculate in explicit 
way an effective boson-boson interaction also when the 
on-site fermion-fcrmion interaction is taken into account. 

The basic idea relies on sharing tasks between the 
fermions and the bosons. In particular, the fermions 
are tuned through a density wave instability estab- 
lishing crystalline order (DLRO), while the bosons 
provide the off-diagonal long range order (ODLRO). 
The interaction between bosons and fermions, and the 
interaction between the fermions, produce an additional 
density modulation also in the bosonic density field, 
hence resulting in a SS phase. To triggering a Density 
Wave (DW) instability in the fermions, the mixed 
boson-fermion system is confined to two dimensions 
and loaded in an optical lattice providing perfect Fermi 
surface nesting at half-filling [24(. 

The Supersolid transition triggered by the fermions 
competes with an instability towards Phase Separation 
in the boson system. Because of the dimensionality 
and the lattice geometry of our system, the presence 
of Van Hove singularities ([25[ and [l6j]) -well studied 
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in BCS superconductivity [26]- strongly enhances the 
tendency towards Phase Separation and produces new 
and interesting features in this transition: an arbitrary 
weak interaction between the bosons and the fermions, 
and between fermions, is sufficient to drive the Phase 
Separation at low temperatures. In the following, we 
investigate the instabilities towards Phase Separation 
and Density Wave formation. We focus on the weak 
coupling limit between the bosons and the fermions 
and between fermions, which excludes a demixing in a 
repulsive Bose-Fermi system along the lines discussed in 

The organization of paper is as follows. In the section 
II, we set the notation, and derive the model Hamilto- 
nian. In the third section, we display the novel result 
consisting in the analytic formula for the effective boson- 
boson interaction as a function of the temperature, and 
discuss the competition between PS and SS; we derive the 
Phase separation and Density- Wave transition tempera- 
tures. In the last section, we write down our conclusions 
and comment about future perspectives of the topic. 



II. THE MODEL HAMILTONIAN 

The microscopic Hamiltonian for interacting spinless 
bosons and not spin-polarized fermionic atoms subject 
to an optical lattice reads H — Hb + Hp + Hi nt (here 
and in the following, a = B,F and a = \A denote the 
atomic species and the spin, respectively) [3, Hf| 



Hi, 



dr 



(1) 



Here, we assume a repulsive interaction between the 
bosons: gsB — a B B /r>iB > 0, where asB is asso- 
ciated the s-wave scattering length. The boson-fermion 
interaction strength qbf = ^ttU a B p Ithr is assumed 
to be the same for both the spin configurations; here, 
clbf is the boson- fcrmion s-wave scattering length and 
to/? = mBfUp I {tub + tif) is the reduced mass. Finally, 
gFF — 2i:h 2 app/mp with a^p the fermion-fermion s- 
wave scattering length. In the following we always as- 
sume both clbf > and app > 0. 



The optical lattice with wave length A provides 
an a = A/2- periodic potential for the bosons and 

fermions with 



Vb.f(x,u) = V b .f( sin 2 — 

a a ' 

[291 ] . where the lattice depth Vp is assumed the 
same for both the spin configurations. In the ex- 
periments, the 2D setup is realized by generating 
an anisotropic three-dimensional optical potential 



sin 



2 ^y, 



V B ,F(x,y) = V b ,f (sin' 



sm 



V, 



7T Z 



73, F 



sm 



with V, 



73, F 



> ^T3.f; 



a a a z 

then the interplane hopping is 



quenched. 

Due to the strong localization around each lattice site 
fi, the annihilation (creation) bosonic and fermionic 
field operators \I/ Q (^J,) may be expanded in terms of 
the Wannier functions w l a (r — r^), with I the Bloch band 
index (H 



n) , 



(2) 



where d\ is the bosonic (&') or fermionic (c l ai ) annihila- 
tion operator acting on a particle at the ith lattice site 
and in the Ith Bloch band. For a strong optical lattice 



Vb,f > E r B F = 2h 7r /A m.B,F the restriction to the 
lowest Bloch band (I = 0) is justified pj. The transla- 
tionally invariant lattice version of Hamiltonian ([T]) reads 



H = -J B blbj-Jp ^2 4.o- g j,o-+ ^ ni(nj-l)+UBF y^(ntWj, -)+E/ r F ,F y^m^m^jj+a y^m^- rhj^) 

<i,j> <i,j> .a i i,cr i i 

(3) 



where we have omitted the band index 1 = 0. The sym- 
bol < i, j > denotes couples of nearest-neighbor lattice 
sites, and 5 the imbalance between spin-up and spin- 
down fermions; hi = b\bi and rhi t „ = c\ a Ci^ a are the num- 
ber operators for bosons and fermions with spin a at the 
ith site, respectively. The boson-boson, boson-fermion, 
and fermion-fermion interaction amplitudes are L^bb = 
.97373 / dr\w B {r)\^, Ubf = gBF J dr\w B {r)\ 2 \wp(r)\ 2 ^ and 
U FF = gFF J dr\w F ,i(r)\ 2 \wp.i(r)\ 2 , respectively [H, H, 




Among the adavantages provided by optical lattices, 
there is the possibility of tuning the Hamiltonian param- 
eters in such a way to realize different interaction regimes. 
Here we focus on weak coupling regime, which takes place 
when \ BF = U 2 BF N a /U B B « 1 (N = 1/2tt 2 J f ) and 
ts = 8Jb/tibUbb >> ^73F ; where ns is the bosonic fill- 
ing factor. 

The parameters involved in the Hamiltonian ^ are re- 
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lated to characteritistic quantities of the optical potential according to 



J 



Jb.f = (4V5F)^, F Vj£exp(-2v^V); 
^f~ =8 ^ 1 + VVf/V b — (Vf) {Vf) > 

A7 



J B 

Uff 
EL 



A7 



(4) 



r 



where 7 = 2a z /X. The hopping and the atom-atom 
interaction amplitudes (0]) are evaluated by extending to 
our case the calculations performed in [l], y, • 



By exploiting the discrete Fourier transform of a,i and 
of its Hermitian conjugate, the Hamiltonian ^ may be 
written in the momenta space 



U2 E [ 



(5) 



r 



where M is the number of lattice sites in x (y) 
direction. We assume the wave vector k be re- 
stricted to the first Brillouin zone: k x G [— it/a, ir/a], 
k y £ [— 7r/a, n/a]. The density number operators are 

n B ,k = Sp^+pV K,n = 'Ep^+f^P 0; the 

bosonic and fermionic dispersion relations e B r and 
— 2Js(2 — cos(fca;a) — cos(fc y a)) and 



£ ff ,fc read e B,k 



-2JF t n{cos(k x a) — cos(k y a)) [28 



III. PHASE SEPARATION VERSUS 
SUPERSOLID 

In this section, we want to gain analytical insight into 
stability of the mixture ground state. As explained in 
to achieve such a goal, the starting point is the 
derivation of an effective Hamiltonian for the bosons 
alone, by keeping in mind that fermions enters the de- 
scription of our system via a modified interaction between 
the bosons themselves. Within linear response theory 
[H, [l9[ the boson density n_e(g) (for the sake of simplic- 
ity, we refer to nB(o)(fn <r(#)) as to the induced pertur- 
bation) drives the fermionic system trough the following 



modulation of the density 

< m a £ >= X<y{UBFn B ,q + U FF < m_ a £ >) 



(6) 



The function \ a , for the a component of the spin, is the 
Lindhard function depending upon the absolute temper- 
ature T and on the wave vector q 

[IE HI: 



1 / , rffc/(e (T ,g)-/(e CT| g +f ) 



Vo 



(7) 



a,k-\-q 



where vq = (27r/a) 2 is the volume of the first Brillouin 
zone; the integration is performed over this region. 
The temperature T enters via the Fermi distribution 



function f(e g g) = !/((! + exp(- 



-)) , with /ip the 



chemical potential of the fermionic atoms. 
We analyze the behavior of the system for temperatures 
well below the the superfluid transition temperature Tkt 
(Kosterlitz-Thouless) . Hence, the mixture is enhanced 
in a sufficiently low temperature regime so that the 
fermionic chemical potential /if my be safely identified 
with the Fermi energy Ep. 

Proceeding form the Hamiltonian ([5]) , we integrate out 
the fermionic freedom degrees following the same path as 
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in [l2|. Within the procedure of tracing out the fermions, 
we treat the spin-up component independently from the 

I 



other, i.e. we perform a mean-field approximation. We 
get the Hamiltonian 



*%Tf = Jp I E -Y'^B^B.-k + u BF{n Bj : < m ? _ £ > +n B i < rh u _ % >) + U FF < >< rh u _ % > ] , (8) 

k 



r 



which describes the effective interaction between the rized in Eq. ([6]) in the Hamiltonian JSJ, we obtain 
bosons of the mixture. By employing the rules summa- 

I 



Hill = E [ U BF{n B ^U BF h B _- k + x^U FF < m u _ % > ) + 

k 

U BF{fi B j:{xiUBFn B _j: + X[Uff < ro T _g > ) + 
UFF(xiU B Fn B £ + X]Uff < rh^ >){xiUBFri B _g + xi^ff < ^ T _fc >)]■ 

I 



(9) 



We exploit the rules ([6]) in Eq. ([9]) in iterative way. Then, 
the Hamiltonian §§§ can be expressed as an expansion in 
series of powers of Ixo-TffI; we have verified that this 
last quantity is much smaller than one. We assume to be 
in a situation in which the spin-up population is equal to 
the spin-down one (5 = 0), and in which the fermions are 
in half-filling configuration, Ep = and = mi = 1/4. 
In such a situation the series associated to the Hamilto- 
nian ((9]) may be summed. We employ the fact that when 



J 



8 = 0, is e-f = = cf and then X] — Xl = Xi with 

^^i^/Mz/W, (10) 



i'O 



£ F,fc € F,k+q 



Together with these two last properties, we use the sym- 
metry e F £ = e F „jr; then, the effective boson-boson in- 
teraction ([5]) reads 



rrint 

eff ~ 2M 2 



E 



U BB + 2 X (T,q)U 2 BF [ 



Uff 



1 - X (T, q)U FF (X{T, q) 2 U FF - 1) 



(11) 



r 



The effective boson-boson interaction U e ff(T) depends 
on the temperature according the formula 



U eff (T) = U BB +2 X (T,^U 2 BF ( T 

X(T,q)U FF s 
(x(T,q)UFF-l) 2h 



1 



x{T,q)U F F 
(12) 



The novelty of the present work relies on Eq. (|T2"j) . Such 
a formula represents a very useful tool for investigating 
in analytical way the instability related to the Van Hove 
singularity and the one associated to the Density Wave. 
In particular, we are interested in calculating the tem- 
peratures of transition to PS and to DW phases. 
To this end, let us focus on fermionic dispersion relation 



e F £ = — 2 Jp[cos(k x a) + cos(fcj,a)]. By analyzing e F g, we 
realize that the Lindhard function (fTUJ) exhibits two loga- 
rithmic singularities. These two singularities give rise to 
instabilities in the system. In particular, the singularity 
at q = induces an instability towards Phase Separation. 
On the other hand, in correspondence of the wave vector 
kuw = (ft/ a, ft/a) -which joints the saddle points (SP) 
ksp = (0, ±7r/a), (±n/a, 0) of the fermionic dispersion 
relation- the fermionic energy vanishes. The wave vector 
kow drives a Density Wave, responsible for the Super- 
solid order. 

Let us focus on the logarithmic Van Hove singularity, and 
analyze the density states. Due to the balance between 
the spin-up and the spin-down populations, we have that 
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A^-(e) = iV| (e) = N(e). Within the tight-binding regime, 
N(e) reads [ij, [3 



w(e , = Wi- (e±<)2 ' 



N , , 16 J F , , 
-In |^|, (13) 



with K[k] the complete elliptic integral of the first kind 

& 

We write the response function (|10p in the energy space 
[T(ij |. and analyze it at <f = and in the limit of zero 
temperature 



X (T^0,0) = 



HJp 



deN(e)s=od e f(e) EF =o 
--fln(^), (14) 



where 8J F is the bandwidth, and the subscripts 5 = 



J 



and Ep = denote the absence of population imbal- 
ances and the half-filled band situation, respectively; 
ci = (2exp(C))/7r w 1.13 with C the Eulero constant 
(« 0.577...). The coupling between the bosons and the 
fermions induces an attraction between the bosons, and 
the effective long distance scattering parameter takes the 
form 



Ueff(T) = U t 



BB 



2 X (T,Q)U 2 BF ( 



1 



1- X (T,0)U FF 
X (T,0)U FF 



1^)-(15) 



(X(T,0)U FF -1) 

The thermodynamic stability of a superfluid condensate 
at low temperatures requires a positive effective inter- 
action, i.e. U e ff(T) > as commented in [31]. The 
condition U e ff(Tps) = defines the critical temperature 
Tps for Phase Separation 



Tps = l6c\J F exp[ 



N Q U FF - (1 



r 



j _ 2N U FF \ \ 



(16) 



J 



Let us focus, now, on the instability in the system triggered by kp>w- By exploiting the symmetry e F ^ + % L 



the energy representation of the response function, in the energy representation, reads 

X(T, k DW ) 



-2e 

■ — (ln(— y— )) • 



Within Bogoliubov theory, the bosonic quasiparticle spectrum is [3l| 



- Jc B (g) 2 + 2n B e B (q)[U eff (T)}. 



(17) 



(18) 



The induced attraction between the bosons reduces the energy of quasiparticle providing a roton minimum at kp)Wi 
which vanishes at critical temperature Tp>w given by 



Tdw = 16ci J F exp 



1 



I 



2Ar 



-N {2 + t B )U FF +X BF 



(19) 



2+t f 



By performing the limit U F p — > in Eq. (|16[) and in Eq. 
(|19|) we are able to reproduce the temperatures signing 
the Phase Separation and Density Wave calculated in 

The fermion-fermion interaction plays a crucial role in 



determining the order characterizing the ground state of 
the mixture. The interaction between the fermions in- 
fluences the behavior of the system via the expression 
of the effective boson-boson interaction, as summarized 
in Eq. (|12p . A non vanishing repulsive fermion-fermion 
interaction acts in such a way to rising a barrier of poten- 
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FIG. 1: On the horizontal axis the fermionic optical depth 
in fermionic recoil energy units; on the vertical axis the tem- 
perature normalized to the fermionic recoil energy. The dot- 
dashed and continuous lines are the PS and DW tempera- 
tures, respectively, when Uff > 0. The curves are plotted for 
(iff = 104.8 ao, with ao the Bohr radius, 7 = 3, ns = 3/2, 
V§ = 20,V b /Vf = Vg/V£ = 7/3. 



tial, which prevents system to achieve Phase Separation. 
This reflects in a mechanism which advantages the Su- 
persolid order, as shown in the plots reported in Fig. 1 
(Uff > 0) and in Fig. 2 (Uff = 0), that are represented 
for a mixture made up of bosonic atoms of 87 Rb and of 
fermionic atoms of 40 Rb in a square optical lattice with 
A equal to 830nm. From these figures, we see that the PS 
and SS lines cross in a certain point, characterized by a 
critical depth. We observe that the PS region of Fig. 1, 
plotted for Uff > 0, is smaller than the same region of 
Fig. 2 obtained with Uff = 0. The difference between 
the magnitudes of the two PS zones may be interpreted 
as a signature of the role of the finite fermion-fermion 
interaction. The amount of the shift of a PS region with 
respect to the other might be used as an indirect measure 
of the fermion-fermion s-wave scattering length. 
The role of the temperature is important in determining 
the kind of order that in the system establishes as well. 
In fact, when the optical lattice depth Vf is greater than 
its critical value, the greater is the temperature the more 
the Phase Separation is dominant in the system; below 
the critical depth, the behavior is reversed. The compar- 
ison between the absence of the on-site fermion-fermion 
interaction and the case in which such an interaction is 
present, is shown in Fig. 3. Here, we observe that when 
the fermions of mixture interact with each other, the SS 
region is wider with respect to Uff = case. Then, we 
may deduce that the Fermi- Fermi interaction acts in such 
a way to make more important the Supersolid order. 



IV. CONCLUSIONS 

In this paper, we have investigated the ground state 
of a Bose-Fermi mixture loaded in a two-dimensional 
periodic geometry, and made up of spinlcss bosons 
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FIG. 2: The dashed and dotted lines are the PS and DW 
temperatures, respectively, when Uff = 0. Here is, clearly, 
(iff = 0, and the other quantities assume the same values as 
in Fig. 1. 



y 




3 3.5 4 4.5 5 5.5 



v F 

FIG. 3: The black continuous line represents the separation 
between SS (the points above the line) and the PS (the points 
below the line) when off = 0; the dotted line is plotted in 
with ciff = 104.8 ao- The other quantities assume the same 
values as in Fig. 1. 



and not spin-polarized fermions. We have analyzed 
the competition between Phase Separation and Super- 
solid orders. This analysis was carried out within the 
framework of the linear response theory. By employing 
an iterative technique, we have obtained an analytical 
expression for the effective boson-boson interaction as 
a function of the absolute temperature of the sample. 
We have performed our study in absence of fermionic 
population imbalances and in presence of boson-fcrmion 
and fermion-fermion repulsive interactions. We have 
studied the phase diagram of the system by stressing the 
role of the fermion-fermion interaction in determining 
the kind of order sustained by the Bose-Fermi mixture 
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Hamiltonian. We have explained the changes that such 
an interaction introduces in the behavior of the system 
with respect to the case in which the fermions do not 
interact on the same lattice site. We have discussed the 
nice interplay between the critical value of the optical 
lattice depth for the fermions and the temperature of the 
sample in establishing the kind of order in the ground 
state of the mixture. 

The study of such a topic opens up a bunch of very 
exciting possibilities. In particular, a very interesting re- 
lationship could be established between these issues and 
the ones related to propagation of the zero sound waves 



both in homogeneous [18[ and in trapped [32J systems of 
ultracold atoms; for these kind of problems, in fact, the 
linear response theory is exploited as well. 

In a forthcoming paper, we wish to accomplish the 
task of analyzing the case of imbalance between the 
two fermionic populations and different sign of boson- 
fcrmion and fermion-fcrmion interactions, by employing 
a suitable modification of the iterative method displayed 
in the present work. 

This work has been partially supported by Fondazione 
CARIPARO. Discussions with S.M. Giampaolo are ac- 
knowledged. 
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